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Q . ABSTRACT 

, In this paper we fit an analytic function to the Bivariate Brightness Distribution 

(BBD) of galaxies. It is a combination of the classical Schechter Function convolved 

I \ with a Gaussian distribution in surface brightness: thus incorporating the luminosity- 

^ . surface brightness correlation as seen in many recent datasets. We fit this function to 

' a recent measurement of the BBD based on 45,000 galaxies from the two-degree field 

\ Galaxy Redshift Survey (Cross et al. 2001). The parameters for the best fit model are 

. (/)* = (0.0206±0.0009)/i3Mpc-3, M;^ -5 log h = (-19.72±0.04) mag, a = -1.05±0.02, 

O ; ^ 0.281 ± 0.007, = (22.45 ± 0.01) mag arcsec'^ and = 0.517 ± 0.006. cf)* , 

I ' M^, and a equate to the conventional Schechter parameters. is the slope of the 

. luminosity-surface brightness correlation, /i* ^ . is the characteristic effective surface 

I brightness at M^, and cr^ is the width of the Gaussian. 

Q^' Using a BBF we explore the impact of the limiting detection isophote on classi- 

' I cal measures of the galaxy luminosity distribution. We demonstrate that if isophotal 

^ . magnitudes are used then errors of AAf^*. ~ 0.62 mags, A0* ~ 26% and Aa ~ 0.04 

\ are likely for ^iiim,bj — 24.0 mag arcsec^^. If Gaussian corrected magnitudes are used 

'. these change to AM^ - 0.38 mags, Acj)* - 11% and Aa < 0.01 for ni„n,b, = 24.0 

^ , mag arcsec" . Hence while the faint-end slope, a, appears fairly robust to surface 

■ brightness issues, both the M* and 0* values are highly dependent. The range 
' over which these parameters were seen to vary is fully consistent with the scat- 
ter in the published values, reproducing the range of observed luminosity densities 
(1.1 < jbj < 2.2 X lO^/iLoMpc"^ see Cross et al. 2001). If total magnitudes are 
recovered then there is no change in the luminosity function within the errors for 
l^iim,bj — 24.0 mag arcsec"^. We conclude that surface brightness selection effects are 
primarily responsible for this variation. After due consideration of these effects, we 
derive a value of jbj = 2.16 x IO^/iLqMpc"^. 
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1 INTRODUCTION 

For the past quarter of a century the luminosity distribu- 
tion of galaxies has been represented by a Schechter func- 
tion (Schechter 1976). This contains three defining parame- 
ters (see Eqn. ^ and from these parameters one can derive 
useful quantities such as the mean local luminosity density. 
Over the past decade many measurements of the Schechter 



parameters have been made (e.g. Efstathiou, Ellis & Peter- 
son, 1988; Loveday et al. 1992; Marzke et al. 1994, 1998; 
Zucca et al., 1997). However when compared, the results 
from these recent surveys show a wide range in the derived 
Schechter parameters, Cross et al. (2001) and hence a wide 
range in the inferred luminosity density. 

The two most likely culprits for this variation are; cos- 
mic variance and/or surface brightness selection bias (e.g. 
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Disney 1976; Phillipps, Davies & Disney 1990). In Cross 
et al. (2001) these two errors were quantified for a sample 
of 45,000 galaxies drawn from the two-degree field galaxy 
redshift survey (2dFGRS). The conclusion was that both of 
these issues are significant and can potentially lead to signif- 
icant underestimates of the local luminosity density (~few 
% for cosmic variance and ~35% for surface brightness bi- 
ases). The surface brightness biases can be separated into 
two distinct parts; a surface brightness dependent photo- 
metric bias and a surface brightness dependent Malmquist 
bias. The biggest effect comes from the former, i.e. underes- 
timating the total magnitudes (25%). The remaining contri- 
bution comes from overestimating the volume over which a 
galaxy can be seen (10%). Cross et al. (2001) advocated the 
construction of a Bivariate Brightness Distribution (BBD) 
to remove the surface brightness selection biases. The BBD 
is the space density of galaxies as a function of the absolute 
magnitude, M, AND the absolute surface brightness, /i. This 
process also revealed a number of noteworthy results in its 
own right: Firstly a clear dearth of giant low-surface bright- 
ness galaxies; secondly a clear luminosity-surface brightness 
correlation (spanning —23 < Mb- < —16); and thirdly the 
general rise in the space-density from giant to dwarf systems 
to the faint limit {Mb- = -16). 

Some earlier attempts have been made to measure the 
BBD (e.g. Choloniewski 1985; van der Kruit 1987; Sodre & 
Lahav 1993). These were for small samples of bright galax- 
ies. The Choloniewski sample contained 248 E/SO galaxies 
with —22 < Mb < —18; the van der Kruit sample con- 
tained 51 galaxies with isophotal diameters > 2' at 26.5 mag 
arcsec"^ (in the photographic Illa-J band) and the Sodre & 
Lahav sample contained 529 galaxies with isophotal diam- 
eters > 1' at 25.6-26.0 mag arcsec"^ in the B-band. These 
are very strict limits, and only very intrinsically bright and 
large galaxies were well sampled. 

Recently three more extensive and independent mea- 
surements of the Bivariate Brightness Distribution have 
been made. Driver (1999) derived the BBD for a volume- 
limited sample drawn from the Hubble Deep Field, de Jong 
& Lacey (2000) studied a homogeneous sample of 1000 late- 
type spirals and Blanton et al. (2001) derived the BBD for 
a sample of 11,275 galaxies from the Sloan Digital Sky Sur- 
vey. These three surveys plus Cross et al. (2001) all confirm 
the existence of the luminosity-surface brightness correla- 
tion, demonstrating that surface brightness selection biases 
are luminosity dependent. The luminosity-surface brightness 
correlation is measured to be Mb^ = 2.4 ±J;| fie — 72.3±32;9 
from (Cross et al. 2001). This equates to fie oc 0.42 ±aie, 
Mt^ -|-30.2±i:l. From Hubble Deep Field data. Driver (1999) 
found a steeper gradient {Mpisow oc 1.5/ie) as did Ferguson 
& BinggeU (1994) in the Virgo cluster {Mb oc 1.4^io). The 
number density of galaxies is a maximum along this line, 
falling away at both higher and lower surface brightnesses. 
Given the luminosity-surface correlation it is hardly surpris- 
ing if surveys with differing selection criterion recover widely 
ranging Schechter parameters. 

A useful next step is to produce an analytical function 
to fit to these derived BBDs. In §2 we show such a function 
and fit it to the derived 2dFGRS BBD. In §3 we discuss 
our fitting procedure and compare it to a similar estima- 
tion made by de Jong & Lacey (2000). In §4 we show how 
the luminosity density can be calculated from the BBF and 



compare it to the published values. Finally we explore the 
scope for error in Schechter function estimators if surface 
brightness selection effects are ignored. 



2 A BIVARIATE BRIGHTNESS FUNCTION 

The luminosity function as traditionally described by the 
Schechter Function (Schechter 1976), is shown below as 
Eqn.|l} 

<P{M)dM = 0.4 In(lO)^* io0.4(M*-M)(.+i) ^_io°-(-*-») 

(1) 

This equation contains three parameters M*, (p* & a 
which describe the "characteristic magnitude", the "nor- 
malisation constant" and the "faint end slope", respec- 
tively. Containing no surface brightness information this 
function provides a good fit to the space density of field 
galaxies albeit over a fairly restricted range of luminosities 
(-22 < Mb- < -16) in the field. However, the Schechter 
Function provides a poorer fit to the luminosity distribution 
in clusters (e.g. Driver et al. 1994, Andreon, Cuillandre & 
Pello 2000) and when the population is subdivided according 
to spectral type (Madgwick et al. 2001). Current estimates 
constrain the three Schechter parameters, for field galaxies, 
to lie in the range: -19.75 < M* < -19.15; 0.013 <<;/>*< 
0.027; —1.22 < a < —0.7 resulting in a pessimistic luminos- 
ity density range of: 1.1 < jb- < 3.2 x 10*/i LoMpc"^. 

Empirically Cross et al. (2001) find that the distribu- 
tion of the galaxy population in surface brightness appears 
symmetrical about a ridge and can therefore be described by 
a Gaussian or possibly a quadratic distribution. The ridge is 
described by /i = /9^M -I- C, the luminosity-surface bright- 
ness correlation. Clearly a BBF needs to relate closely to 
the Schechter function in luminosity. Hence by multiplying 
the classical Schechter function with a Gaussian in surface 
brightness we can construct the following BBF: 

1 , fie - fl*e - 13^ {M - M,) ^2, 

''''P^ 2^ ' J (2) 

where /J^ is the gradient of the luminosity-surface bright- 
ness correlation and o-p is the dispersion in the surface 
brightness. This function is identical to that presented by 
Choloniewski (1985), and derived by Dalcanton, Spergel & 
Summers (1997) and de Jong & Lacey (1999a,b 2000) using 
the Fall & Efstathiou (1980) disk-galaxy formation model. 
Note that the new term contains a normalisation coefficient 
■^y= — ensuring that </!>*, M* and a are identical to the tra- 
ditional Schechter parameters. 



3 FITTING THE BBF 

We choose to fit the BBF to the data shown in Cross et 
al. (2001); this is the largest available data set. This BBD 
was derived from a subset of 45,000 galaxies from the two- 
degree field galaxy redshift survey (see Cross et al. 2001 for 
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Figure 1. The thick Uncs show the BBF computed for the best fit 
parameters. The thin lines depict the 2dFGRS BBD (q.v. Cross 
et al. 2001 Fig. 10b). The contours are at 1.0 X 10"'', 1.0 X 10"^, 
2.5 X 10^3^ 5 X 10-3, 7.5 x 10"^, 1.0 x 10"^, 1.25 x 10^^, 

1.5 x 10-3, 1.75 X 10-2, 2.0 X IQ-^ galaxies Mpc-^ mag-^ (mag 
arcsec-2). The shaded area shows the selection boundary - see 
Cross et al. (2001) for details. 

details) .and is shown on Fig. ^ (thin contours). Also shown 
in Fig. y is the selection boundary derived from visibility 
theory (shaded region, see Appendix B of Cross et al. 2001). 
In the shaded region insufficient volume is surveyed to make 
any meaningful statement of the space-densities. 

The BBF can be fitted to the BBD by minimising the 
of the model compared to the data. The BBF is a non-linear 
six parameter equation and to find the minimum, we use the 
Levenberg-Marquardt Method (see Press et al. 1986). The 
data provided (see Cross et al. 2001) is binned as in Table 
C2 of Cross et al. (2001) using those bins whose values are 
binned on a minimum of 25 galaxies. 

The best fit parameters we derive are: </>* = (0.0206 ± 
0.0009)/i3Mpc-3, M* - Slog ft = (-19.72 ± 0.04) mag, 
a = -1.05 ± 0.02, = 0.281 ± 0.007, = (21.90 ± 0.01) 
mag arcsec-^ and — 0.517 ± 0.006. All errors are la 
errors. Fig. |l] shows the BBF for these parameters (thick 
lines) overlaid on the data (thin contour lines) . The errors in 
the parameters were found using a Monte-Carlo simulation, 
that is the observed distribution was randomised within the 
quoted Icr errors and the BBF fit re-derived. The final BBF 
fit yields a value of 164, for v = 49, where v is the no. 
of data points - no. of parameters. This gives a likelihood of 

2.6 X 10-^*. 

Hence, although the BBF appears to describe the BBD, 
the fit is poor. It is important to understand where the dif- 
ferences are occurring. From Fig. ^ we see that the model 
fits the data well brightwards of M = —18, and less well 
in fainter bins. The errors become comparable to the space 
density faintwards of M = —16, so the main error in the fit 



occurs in the range: — 16 > Af > —18. The data (Fig. |l|) 
show an upturn towards the faint end in this range whereas 
the Schechter function gradually flattens towards the faint 
end. Thus it is the Schechter function part that does not de- 
scribe the data well. Note that the BBF provides Schechter 
parameters comparable to the range from previous surveys. 

The model fits the data well in the surface brightness di- 
rection, implying that a Gaussian distribution is a good de- 
scription of the space density as a function of surface bright- 
ness, for a constant absolute magnitude. 

3.1 Comparison with de Jong & Lacey and 
Blanton 

Table |^ compares our BBF with the de Jong & Lacey (2000) 
BBF which was determined for Sb-Sdm galaxies only. As de 
Jong & Lacey use 95% confidence intervals for their errors, 
we have quoted 2a" errors rather than Icr errors for our val- 
ues. We converted their half-light radii parameters to our 
effective surface brightness parameters. De Jong & Lacey 
fit disks and exponential bulges to their data, taking into 
account inclination and internal extinction. In Appendix A, 
we estimate the uncertainty in our results due to not tak- 
ing this into account and find that the error is ~ 0.1 mag 
in Af* and ~ 0.55 mag arcsec-"^ in /Xe- We note that the 
de Jong & Lacey (2000) parameters are their total galaxy 
parameters, not their disk-only galaxy parameters. In each 
case, these have been converted to fe^-band magnitudes us- 
ing B ~ I = 1.7 mag from de Jong & Lacey (2000), and 
bj = B - 0.28(B - V) (Maddox, Efstathiou & Sutherland 
1990), using a value of [B — V) = 0.5 for a late type spiral 
(Coleman, Wu & Weedman, 1980, Driver et al. 1994). Fi- 
nally we convert from Hq = 65 km s-^Mpc-^ to Ho — 100 
km s-^Mpc-^. In addition, the de Jong & Lacey sample 
has tighter selection criteria and more accurate CCD pho- 
tometry (±0.05 mag compared to ±0.2 mag), but it only 
includes late-type galaxies and has a redshift completeness 
of 80% compared to > 90% for the 2dFGRS. In spite of this, 
we find a similar spread in ^ (cr^ = 0.52 q.v. 0.61) and we 
find that the a values of the two surveys are equal within the 
errors. The 2dFGRS has a brighter M* by 0.05 mag, and 
a brighter /i* by 0.85 mag arcsec-"^. Taking into account 
the effects of bulges and inclination as mentioned above, the 
2dFGRS distribution has become 0.15 mag brighter than the 
de Jong & Lacey distribution and has a brighter /x* by 0.3 
mag arcsec-"^. Considering that late-type galaxies tend to 
be fainter and lower surface brightness, these results appear 
fully consistent. 

Our value of can be converted to a luminosity- 
scale size gradient = —0.360 ± 0.004. Although this dif- 
fers from the de Jong & Lacey (2000) value, interestingly, 
it agrees more closely with their theoretical prediction of 
/3re = — ^ (see de Jong & Lacey 2000). One possible reason 
for the variation in may be a correlation between colour 
and absolute magnitude. Blanton et al. (2001) find a strong 
correlation between (g* — r*) colour and Mr* : brighter galax- 
ies are redder, fainter galaxies are bluer. Making estimates of 
the colour-magnitude correlation from Fig. 13 of Blanton et 
al. (2001), we find that Af^. = -8.75±| (5* - r*) - 14.88±2. 
Using a mean bj — r* = 1.1 (calculated from Fukugita, 
Shimasaku & Ichikawa 1995 and Maddox, Efstathiou & 
Sutherland 1990), and assuming that the additional colour 



© 0000 RAS, MNRAS 000, 000-000 



4 Cross & Driver 



term A(&j — r*) = /S.{g* — r*) we calculate that the ex- 
pected Pfi,r* ~ 0.36 ± 0.23. The value estimated from Fig. 
10 of Blanton et al. (2001) is P^,r' = 0.50±o:?. Thus the 
r* band luminosity-surface brightness correlation appears 
steeper than the bj band luminosity-surface brightness cor- 
relation. A similar colour-magnitude correlation in [bj — I) 
could explain the discrepancy between our result and de 
Jong & Lacey (2000) result. 

The general good overall agreement between these sub- 
stantially different surveys is an important vindication of 
both results. Cross et al. (2001) has extended the de Jong 
& Lacey conclusions to the full range of galaxy types with 
M < —16. However, the different values obtained for the lu- 
minosity surface brightness correlation may reflect a colour 
or morphologically dependent luminosity- surface brightness 
correlation. Blanton et al. (2001) seem to have found simi- 
lar results, but have not fitted a function or tabulated their 
results. 



4 CALCULATING THE LUMINOSITY 
DENSITY 

As for the Schechter function it is trivial to calculate the 
luminosity density, j, by integrating the product of the BBF 
and the luminosity over the complete range of surface bright- 
ness and absolute magnitude. 

oc poc 



j = I I L(M)(t>{M,fi)dMd^ 
= 0* r r(Q + 2) = io-oa(m'^Mq) ^ 2) (3) 

The solution is the same as the solution to the inte- 
gral obtained from the Schechter function. When calculated 
using the best fit parameters, the value of the luminosity 
density, jt- = (2.16 ± 0.14) x 10® /iL© Mpc"^. 

In Blanton et al. (2001), the Sloan team get a 40% 
higher value for the luminosity density in the bj filter than 
the 2dFGRS team. Does this mean that 2dFGRS is miss- 
ing some galaxies, or at least underestimating their fluxes? 
For a start the values of M* are consistent, suggesting that 
both surveys are correcting magnitudes properly. However 
the measurement of 0* is over 30% higher in Blanton et 
al. (2001). A more recent paper (Yasuda et al. 2001) re- 
vises the SDSS luminosity density of Blanton et al., to 
jbj = 2.43 ± 0.21 X 10"^/i^Mpc"^. This revision is based on 
a fit to the gala:xy number counts, suggesting that the Blan- 
ton et al. region was overdense by 30%. The revised <f>* value 
{(f)* = 2.05 ± O.l2to!^l X 10"^/i^Mpc"^) is now consistent 
with our measurement of 0* = 2.06 ± 0.09 x 10"^/i^Mpc"^, 
in the bj band. 

Given this revised value of </>*, the Blanton result is 
still 10% higher, but Blanton used a colour term bj — B — 
0.35(5 — V), whereas the correct colour term for the APM, 
tested using EIS data is bj = B - 0.28(5 - V) (Peacock, 
private communication). When these two factors are taken 
into account, the luminosity densities are entirely consistent. 



Cross et al. (2001) used a value of ^0,6 = 5.48, rather than 



the correct value of L 



0,6,- 



5.30. The values for the luminosity 



density in this paper use the correct value. 



As demonstrated in Cross et al. (2001) the peak of 
the luminosity density lies well inside the selection bound- 
aries. When the function is integrated over the range —24 < 
M < —15.5, 20.1 < < 24.1, the value obtained is 
jb ~ 2.14 X lO^h LqMpc~^ as compared to the summed 
BBD which gives: jbj. = (2.11±0.20) x 10*/i LgMpc"^. These 
are the approximate selection boundaries, so the correspon- 
dence is excellent. Unless the distribution shows an upturn 
outside the selection boundary the 2dFGRS data have un- 
covered over 98% of the local B-band luminosity density. 



5 EXPLORING SURFACE BRIGHTNESS 
SELECTION EFFECTS 

In Cross et al. (2001), Fig. 1 showed the variation in the 
LF as measured from a number of recent redshift surveys. It 
was postulated that the large variation was due to surface 
brightness selection effects. In this section the variation of 
the luminosity function with the limiting detection isophote 
is explored and compared to the range of published luminos- 
ity functions. Throughout, an exponential surface brightness 
profile and an Einstein-de Sitter cosmology are assumed. 

To calculate a derived luminosity function, we start with 
our fit to the BBF. We take into account the overestimate 
in /i* by 0.55 mag arcsec"^, derived in Appendix A, and use 
a value of fil = 22.45 mag arcsec"^ [Note that in Cross et 
al. (2001) we used a less sophisticated method to estimate 
the offset caused by the bulge and arrived at a figure of 0.55 
mag arcsec"^]. 

We multiply our updated BBF by a visibility volume 
(Cross et al. 2001) to construct an apparent observed num- 
ber in M and /i (see Fig. The parameters adopted 
to derive the visibility surface are: m faint = 20.0 mag, 

mbright = 14.0 mag, dmin = 2.0", dmax = 250.0", Zmax = 

0.5, Zmin ~ 0. The solid angle used was 300n°. These pa- 
rameters are typical of the observed ranges for the most 
recent surveys. The only parameter allowed to vary is the 
detection isophote which took the values 26, 25, 24, 23.5 and 
23 mag arcsec"^. Fig. ^ shows the complex Malmquist bias 
for fiiim ~ oo, 26, 24, 23 mag arcsec"^ (top to bottom). The 
shaded region shows the approximate location of the galaxy 
population. The lines are contours of constant volume for 
that fiiim- As the limiting isophote becomes brighter, the 
surface brightness dependency of the Malmquist bias in- 
creases. 

Given the observed distribution, each galaxy within 
each bin is then randomly assigned a volume out to the maxi- 
mum derived from visibility theory. This volume is converted 
to a redshift assuming an Einstein-de Sitter cosmology and a 
standard k-correction. We assume that the number density 
does not vary within the bin as a function of redshift, i.e. 
there is no evolution and no clustering. The exact absolute 
magnitude and surface brightness value is then randomly 
assigned within each bin, (this assumes that the number 
density does not vary within the bin as a function of M 
or fie)- Around the Al* point particularly, this assumption 
fails, but a Monte Carlo simulation done at higher resolu- 
tion finds no significant differences. The numbers in Table 2, 
and the plots in in Fig ^ ^ were produced from this Monte 
Carlo simulation. 

The net result is a magnitude-limited sample with ob- 
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Figure 2. This figure shows 3 contour plots. The thin Unes show 
the contours of the visibility surface, in Mpc^bin"'^. Each line is a 
decade apart. The five lines shown range from 10^ to 10^, from the 
bright end to the faint end respectively. This visibility surface has 
a detection threshold of 26.0 mag arcsec"-^ and the limits are in 
isophotal magnitudes (see §5). The medium thickness lines show 
the BBF from Fig. Q offset in surface brightness by 0.55 mag 
arcsec~^. The thick lines show the number of galaxies detected 
in each bin. The contours levels are at 1.0, 3.2, 10, 32, 100, 320, 
1000, 3200 and 10000 galaxies mag"! (mag arcsec-^)-!. 

jects randomly distributed within their allowed volume. We 
now calculate each galaxy's isophotal magnitude, a Gaus- 
sianQ corrected magnitude and their total magnitude and 
sum the final number distribution according to absolute 
magnitude. This is plotted in Fig ^ for total, corrected and 
isophotal absolute magnitudes. 

We then reconstruct the luminosity function using a 
l/Vjvfaa: prescription (as our simulations contain no cluster- 
ing this should be an optimal estimator). Fig. ^ shows the 
recovered luminosity functions. The LFs of Fig. ^ demon- 
strate the impact of surface brightness selection as they are 
all drawn from the same BBF; the only difference is the lim- 
iting isophote and the choice of magnitude measurement. 
The range of published values is shown as the shaded area 

t One simple and popular method to correct for light lost is the 
Gaussian correction employed by Maddox, Efstathiou & Suther- 
land (1990) on the APM, and used as a correction in the Source 
Extractor code (Bertin & Arnouts 1996). It works by fitting a 
Gaussian with central surface brightness, fio, and standard devi- 
ation, (T, to the light profile of the galaxy, such that the isophotal 
radius of the Gaussian matches the isophotal radius of the galaxy 
and the isophotal magnitude of the Gaussian is equivalent to the 
isophotal magnitude of the galaxy. The Gaussian corrected mag- 
nitude is then the total fiux under the Gaussian. This works well 
for compact objects such as stars and small angular scale size 
galaxies where the seeing dominates the profile. 



yU„^=™mags/sq arcsec /i|i^=24 mags/sq arcsec 
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Figure 3. This plot shows contours of equal volume for different 
selection functions. The top, left plot shows the selection func- 
tion for an infinite threshold. The selection lines have no sur- 
face brightness dependence. The other plots show isophotal (thin 
lines) and total (thick lines) selection functions, for the detection 
thresholds 26, 24, & 23 mag arcsec"^. As one goes to brighter 
thresholds, the volume becomes a stronger function of surface 
brightness. Thus the mean volume at any M decreases, and a vol- 
ume correction or estimator which is only magnitude dependent 
becomes more biased. The shaded parallelograms represent the 
number density distribution, in terms of the luminosity-surface 
brightness correlation, and the width of the surface brightness dis- 
tribution. Where the contour lines are vertical, an estimator with 
magnitude dependence only is unbiased, elsewhere it is biased. 
Where the contour line crosses the distribution, the luminosity 
function can be recovered using an estimator which takes into ac- 
count surface brightness. Where the distribution is missed, there 
is input catalogue incompleteness, and so there is not enough 
information to recover the luminosity function. When total mag- 
nitudes are used: at 26 mag arcsec"-^, there is no surface bright- 
ness dependency within the shaded region for A/ < —12, so a 
magnitude-only will give an unbiased luminosity function; at 24 
mag arcsec"^, it will be unbiased for M < —14.3, but an estima- 
tor with surface brightness built in, such as the method in Cross 
et al. (2001), will be unbiased M < —12; at 23 mag arcsec"^, the 
magnitude only estimator will be biased at all magnitudes, but a 
surface brightness estimator will recover the luminosity function 
for M < —16. Using isophotal magnitudes, all luminosity function 
estimators will be biased if they are only magnitude dependent. 



(excluding the LCRS Lin et al. 1996). Also shown is the 
limit solution, for our model BBF. The left panel assumes 
isophotal magnitudes are measured, the central panel as- 
sumes Gaussian corrected magnitudes were used and the 
right panel assumes some procedure has been implemented 
to recover the total magnitudes. The results are also tabu- 
lated in Table ||. 
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Figure 4. This is a plot of the number of galaxies detected as a 
function of absolute magnitude, for a limiting isophote of 24.0 mag 
arcsec"'^. The solid line shows the plot for total magnitude, the 
dashed line shows the corresponding corrected magnitudes and 
the the dotted line shows the corresponding isophotal magnitude. 
The difference between each line is a combination of an offset in 
magnitude and fewer galaxies being detected going from total to 
isophotal magnitudes. The peak of the total magnitudes is ~ 0.3 
mag brighter than the peak in the corrected magnitudes which is 
~ 0.35 mag brighter than the peak in the isophotal magnitudes. 

Isophotal magnitudes 

If isophotal magnitudes are adopted and the surface bright- 
ness limit is bright, the luminosities of galaxies are severely 
underestimated. Thus both the number density and the M* 
value are severely underestimated (see Table H) • The varia- 
tion in (j)* is upto 50% and in M* upto 1.0 mags. This tallies 
well with the range of Schechter values recovered (see §2) 
over the range tested (23 < < 26). To some extent is it 
surprising that 0* is not more drastically effected; this is be- 
cause the observed distribution of galaxies is skewed towards 
the faint end, see Fig. ^. As a simple 1/Vmax correction or 
maximum likelihood estimator based on the isophotal mag- 
nitudes alone does not take into account surface brightness 
issues, especially light loss, a smaller volume is calculated 
than for total magnitudes, leading to an overestimate of the 
number density, see Fig |^. This is tempered by a lower num- 
ber density at brighter absolute magnitudes. 

Perhaps most surprising is the robustness of the faint 
end slope whose value is recovered correctly regardless of the 
isophote. 

Corrected magnitudes 

Most surveys attempt to correct their isophotal magnitudes 
to total magnitudes. We used a Gaussian correction as de- 
scribed above. Fig. and Table S demonstrate that corrected 



magnitudes recover 63% of the luminosity density at 24 mag 
arcsec"^ compared to the 43% that isophotal magnitudes re- 
cover and the 98% that total magnitudes recover. As with 
isophotal magnitudes, corrected magnitudes give a luminos- 
ity function biased at all values of M, although the bias has 
been significantly reduced. 



Total magnitudes 

If some method is employed to correct the galaxies to total 
magnitudes (e.g. Kron magnitudes or Petrosian magnitudes) 
we find that the parameters are very robust for jj,iim > 24 
mag arcsec"'^. However, at )j,iim = 23 mag arcsec"'^ the num- 
ber density is underestimated throughout the distribution. 
Fig ^ illustrates why this occurs. 

The volume has almost no surface brightness depen- 
dency for the thresholds 24 < /iHm < 26 provided M < — 14, 
but it has significant surface brightness dependency for 
-22 < M < -14 at fj,iim = 23. The bright absolute 
magnitudes are affected because of cosmological dimming 
and the K-correction. Galaxies at the maximum redshift 
of z = 0.5 will have an apparent surface brightness 3 mag 
arcsec"^ fainter than their intrinsic surface brightness (1.7 
mag arcsec"^ due to cosmological dimming and 1.3 mag 
arcsec"^ due to the K-correction, where K{z) = 2.5z). Even 
galaxies at z = 0.25 will be fainter by 1.6 mag arcsec"^. 
Thus a galaxy with central surface brightness of 21.5 mag 
arcsec"^ {He = 22.6 mag arcsec"^) and z = 0.25 will not be 
detected with a threshold of 23 mag arcsec"^. 

Recovering total magnitudes beforehand will give good 
estimators for the luminosity function, provided that signifi- 
cant numbers of galaxies are not missing. This is a particular 
problem if your maximum redshift is very high. However, if 
galaxies are missing because of cosmological effects, rather 
than being too intrinsically dim even at z; = 0, the number 
density can be recovered using a surface brightness depen- 
dent volume correction as discussed in Cross et al. (2001). 

Overall Effects 

Overall the variations recovered in M* , 0* and a between 
simulated surveys with limits of 24 < ^ < 26 (i.e., com- 
parable to existing surveys) are —19.74 < Af* < —19.08, 
0.020 <(!>*< 0.017 and -1.07 < a < -1.05. Fig. | shows 
that the observed variation in Schechter function parame- 
ters has been recovered at the bright end for 24 < ^um < 26 
when either isophotal or Gaussian corrected magnitudes are 
used. However, Fig. ^ demonstrates that the observed vari- 
ation in the faint end slope has not been recovered and is 
testament to the fact that surface brightness selection ef- 
fects do not reproduce all the variation seen in the faint-end 
slope. This supports the suggestion in Cross et al. (2001) 
that the faint-end slope depends more critically on the clus- 
tering correction than surface brightness issues. 

Blanton et al. (2001) find a 0.08 change in the faint 
end slope going from ^e,r* ~ 23.5 to /ie,r* = 24.5, with 
no change at the bright end. Using a {bj — r*) = 1.1, this 
leads to a 6j-band isophote of num ~ 23.5 for fie,r* = 23.5. 
The change in the faint end can be compared to the changes 
seen in Table 2, for total magnitudes. The faint end slope, 
a ~ —1.037 for = 23.5 and a ~ —1.058 for /iHm ~ 24.5. 



© 0000 RAS, MNRAS 000, 000-000 



The BBF 7 




Figure 5. The variation of LFs with for isophotal (top), corrected (middle) and total (bottom) magnitudes, /lum varies from 26 

mag arcsec"^ to 23 mag arcsec"^. The shaded region shows the variation of recent surveys from Sloan (with the new ip*) to APM. The 
LCRS is not shown as it has an additional surface brightness constraint that significantly reduces the faint end. 



This gives a 0.021 change over a similar interval, lower than 
the Sloan result. However, Sloan has a red selected sample, 
which gives a steeper luminosity-surface brightness correla- 
tion (see §3.1). This could account for a greater change in 
a. 

In Cross et al. (2001) we take an isophotally selected 
sample and apply corrections in surface brightness as well 
as absolute magnitude. These corrections imply that we will 
not underestimate M* or <j>* . SDSS calculated Petrosian 
magnitudes before selecting their sample. Petrosian mag- 
nitudes are aperture magnitudes and therefore do not show 
such a pronounced variation with redshift as isophotal mag- 
nitudes. Thus the sample is selected from pseudo-total lim- 
its. 

For deep isophotes, the volume correction has virtually 
no surface brightness dependence, as shown in Fig ^ so the 
number density can be calculated at the bright end trivially 
and is only underestimated at the faint end where some 
galaxies have too a low surface brightness to get into the 
sample. Using either of the techniques outlined in the previ- 
ous paragraph (Cross et al. 2001, Blanton et al. 2001) should 
give accurate values of M^_. — 5 log h and (p* , —19.75 ± 0.05 
mag, (2.02 ±0.02) x 10"^/i^Mpc"^ for 2dFGRS (Cross et al. 
2001) and -19.70 ±0.04 mag, (2.05 ±0.12) x lO'^^^Mpc"^ 
for SDSS (Blanton et al. 2001, Yasuda et al. 2001). The 
caveat is that no correction to total magnitudes is perfect, 
and the corrected magnitudes will tend to have some surface 
brightness dependency as is demonstrated by the Gaussian 
corrected luminosity function. Even when isophotal magni- 
tudes have been corrected to pseudo-total magnitudes it is 



better to use a 1/Vmax or maximum likelihood estimator 
which is a function of both M and fi. 

While a good understanding of visibility theory will ac- 
count for galaxies within the surface brightness limits, galax- 
ies with very low central surface brightness, but bright total 
magnitudes can be missed. These are one source of mis- 
matches between the estimates of a in different surveys. 
Another reason for different estimates is inhomogeneities in 
the space density of galaxies. Surveys looking at different 
parts of the sky will encounter variations in the Large Scale 
Structure. Dwarf galaxies are seen over a smaller volume, 
so they can have large clustering corrections. Differences in 
clustering corrections between different surveys will tend to 
bias Q rather than (p* . 



6 CONCLUSIONS 

We have presented a fitting function for the Bivariate Bright- 
ness Distribution. This takes a similar form to a Schechter 
function in luminosity coupled with a Gaussian distribution 
in surface brightness. The Bivariate Brightness Function was 
fitted to the recent results from Cross et al. (2001) who con- 
structed a BBD for a sample of 45,000 galaxies from the 
2dFGRS. The BBF fits the data well at the bright end, but 
poorly at the faint end. 

We compare the parameters of the fit to the de Jong 
& Lacey (2000) results for late-type spirals. While our re- 
sults broadly agree there are differences which may provide 
clues toward understanding formation and evolution tracks 
of different galaxy types. 
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The BBF can be integrated to yield a total luminosity 
density of ji^ = 2.16 ± 0.14 x 10*/i LoMpc"^. This agrees 
with the luminosity density calculated from the data and 
demonstrates that unless the BBD shows sub-structure or 
a dramatic upturn beyond the selection boundaries the ma- 
jority of the luminosity density in the local Universe has now 
been detected. 

This paper has dealt with the biases that occur if you 
start with an isophotally selected sample and do not ap- 
ply any corrections, or apply a light loss correction without 
properly considering the volume correction. In both cases 
the Schechter parameters can be biased, and the luminosity 
density underestimated. 

Using a BBF we explore the impact of the limiting de- 
tection isophote on classical measures of the Schechter func- 
tion. We demonstrate that if isophotal magnitudes are used 
then errors of AM^^. ~ 0.62 mags, A<j)* ~ 26% and Aa ~ 
0.04 are likely at ^.um.bj = 24.0 mag arcsec"^. If Gaussian 
corrected magnitudes are used these change to AMj . ~ 0.38 
mags, A0* ~ 11% and Aa < 0.01 are likely at iium,bj ~ 24.0 
mag arcsec"'^. If total magnitudes can be recovered then the 
observed luminosity function will be correct within the er- 
rors provided fiiim,b- > 24.0 mag arcsec"'^. Hence while the 
faint-end slope, a, appears fairly robust to surface brightness 
issues both the M* point and 0* are highly dependent. The 
range over which these parameters vary is fully consistent 
with the scatter in the published values which come from 
a variety of surveys with differing selection criterion. These 
parameters produce a range in the luminosity density, jt^ , 
of 0.9 < jbj < 2.2 X lO^h LqMpc~^ again agreeing well with 
the range of published values (1.1 — 2.2 x 10*/i LoMpc"^ see 
Cross et al. 2001). 

When selection effects are taken into account properly, 
the luminosity functions of recent surveys agree very well. 
The 2dFGRS and SDSS luminosity functions give the same 
M* and 0* values, but disagree on the values of a. The 
differences in a are likely to be due to a combination of input 
catalogue incompleteness, filter used and different clustering 
corrections. 

§5 suggests that if a deep isophote is used and total mag- 
nitudes are recovered, then traditional magnitude dependent 
estimators can be used. However, future work at higher red- 
shifts or low redshift work on extremely faint galaxies will 
include many galaxies close to the limits of the detection 
threshold, where surface brightness dependencies are strong. 
The BBD provides a framework that allows one to determine 
whether a bias is present (see Fig. ^) and to correct for it. In 
addition, it produces parameters such as l3fj_ and which 
could be used to place constraints on galaxy formation and 
evolution models. 

Our conclusion is that after a quarter of a century we 
need to upgrade our representation of the space-density of 
galaxies to now include surface brightness. The BBD and 
BBF provide an excellent starting point and should lead 
to more reliable and consistent measurements of the local 
luminosity density as well as providing new constraints on 
galaxy formation models. 



Table 1. Bivariate Brightness Function Parameters 



Parameter 


2dFGRS 


dj & L (2000) 


— Slog h 1 mag 


-19.72 ±0.08 


-19.67 ±0.17 


ct 


-1.05 ±0.04 


-0.93 ±0.10 




0.281 ± 0.014 


0.494 ±0.04 


/I* / mag arcscc"^ 


21.90 ±0.02 


22.82 ±0.19 




(22.45 ±0.02) 






0.517 ±0.012 


0.61 ±0.04 
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APPENDIX A: TESTING THE FACE-ON 
PARAMETERS 

In Appendix A of Cross et al. (2001), we calculated the 
difference between our corrected magnitudes and total mag- 
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Table 2. Table of Schechter parameters for Figure 





Magnitude 


M* - 5 log h 

Oj 


<l,*/h^Mpc-^ 






i6^./lO»hL0Mpc-3 


26.0 


Isophotal 


-19.54 ±0.02 


(1.95 ±0.04) 


xlO" 


2 


-1.057 ±0.01 


1.74 ±0.05 


25.0 


Isophotal 


-19.36 ±0.02 


l^i.ou It U.U4 


X lU 


2 


-1.058 ±0.01 


1.40 ± 0.05 


24.0 


Isophotal 


-19.10 ±0.02 


(^l.oz it U.U4 


X iU 


2 


-1.090 ±0.01 


0.928 ± 0.05 


ZO.O 


Isophotal 


— io.oD It U.Uz 


(l.oD it U.U4 


X lU 


2 


— i.Uoo it U.Ui 


U.DDZ It U.UO 


zo.U 


Isophotal 


1 Q Rft -1- n no 
— io.OD it U.Uz 


{l.VKJ ± U.U4 


X lU 


2 


— i.lzi it U.UI 


U.41 / it U.Uo 


26.0 


Corrected 


-19.64 ±0.02 


(2.00 ± 0.04) 


xlO" 


2 


-1.056 ±0.01 


1.96 ±0.05 


25.0 


Corrected 


-19.54 ±0.02 


(1.96 ±0.04) 


xlO- 


2 


-1.055 ±0.01 


1.75 ±0.05 


24.0 


Corrected 


-19.34 ± 0.02 


(1.83 ±0.04) 


xlO" 




-1.051 ± 0.01 


1.36 ± 0.05 


23.5 


Corrected 


-19.22 ±0.02 


(1.62 ±0.04) 


xlO" 


2 


-1.062 ±0.01 


1.08 ±0.05 


23.0 


Corrected 


-19.00 ±0.02 


(1.43 ±0.04) 


xlO" 


2 


-1.021 ±0.01 


0.760 ±0.05 


26.0 


Total 


-19.74 ±0.02 


(2.00 ±0.04) 


xlO" 


2 


-1.066 ±0.01 


2.16 ±0.05 


25.0 


Total 


-19.74 ±0.02 


(2.00 ±0.04) 


xlO" 


2 


-1.065 ±0.01 


2.16 ±0.05 


24.0 


Total 


-19.72 ±0.02 


(2.00 ±0.04) 


xlO" 


2 


-1.058 ±0.01 


2.11 ±0.05 


23.5 


Total 


-19.66 ±0.02 


(1.81 ±0.04) 


xlO" 


2 


-1.037 ±0.01 


1.78 ±0.05 


23.0 


Total 


-19.62 ±0.02 


(1.08 ±0.04) 


xlO" 


2 


-0.975 ±0.01 


0.989 ±0.05 



nitudes compared to isophotal magnitudes and total magni- 
tudes for 8 galaxy types. We have extended this by taking 
seeing into account and looking at the errors in effective sur- 
face brightness too. However, the galaxies are still face on. 
Fig Al shows the errors in magnitudes for isophotal mag- 
nitudes (triangles) and our corrected magnitudes (squares) 
compared to total magnitudes, for the 2dFGRS data — 
^^Um = 24.67 and seeing of 2". Each is shown for 3 different 
redshifts, 0.05, 0.1 and 0.15. The corrected magnitudes are 
quite good with a typical offset of 0.1 mag and the greatest 
offset is 0.2 mag for an elliptical galaxy with a redshift of 0.1 
or 0.15 and a, z = 0.15 LSBG, with a bulge-to-ratio of 0.12. 
The isophotal magnitudes are much worse, with a typical 
offset of 0.4 mag and the greatest offset is 2.1 mag for the 
z=0.15 LSBG. There is only a weak redshift dependence for 
the corrected magnitudes. 

Fig shows the bias in effective surface brightness for 
our correction, before and after seeing. There is up to 1.75 
mag arcsec"^ offset before seeing is taken into account and 
up to 0.7 mag arcsec"^ after seeing is taken into account. 
The greatest offsets occur in galaxies with large bulges, such 
as Sa/SOs and ellipticals. The average offset will be around 
0.4 mag arcsec"^. 

However, there is a significant bias in the effective 
surface brightnesses. Ellipticals have their effective surface 
brightnesses underestimated by ~ 0.5 mag arcsec"^; SO and 
Sa galaxies have their effective surface brightnesses overes- 
timated by up to 0.5 mag arcsec"^. Late type spirals and 
irregulars will only have a negligible effect. 

However late types will be more effected by inclination. 
To model this, we have assumed that a disk galaxy is opti- 
cally thin and has no internal extinction. 

A galaxy of area Aiso, was assumed to have an isophotal 
radius riso, given by: 




0.3 0.4 0.5 0.6 0.7 
Bulge-to-Total Ratio 

Figure Al. This plot shows the bias in corrected magnitudes 
and isophotal magnitudes as a function of bulge-to-total ratio for 
three different redshifts, z = 0.05, 0.10, 0.15. The squares show the 
offset between exponentially corrected magnitudes (Cross et al. 
2001) and total magnitudes for an elliptical galaxy, B/T = 0.999, 
an SO galaxy, B/T = 0.65, Sa-Sd galaxies, B/T = 0.5 - 0.1, an 
irregular galaxy, B/T = 0.001, and a LSBG with, fio = 23 mag 
arcsec"^ and B/T = 0.12. The offset is low generally around 0.1 
mag, with a weak dependence on redshift. The triangles show the 
difference between isophotal magnitudes and total magnitudes as 
a function of bulge-to-total ratio for the same galaxies. There is a 
much greater difference, which is a stronger function of redshift. 



(Al) 



The isophotal magnitude and radius were used to calcu- 
late the total magnitude and the effective surface brightness 
He as described in Cross et al. 2001. 



However a disk galaxy, inclined at an angle, i with major 
axis a and minor axis b has an area 



A, 



TT ah 



(A2) 



© 0000 RAS, MNRAS 000, 000-000 



10 Cross & Driver 




0.3 0.4 0.5 0.6 0.7 
Bulge-to-Total Ratio 

Figure A2. This plot shows the bias in effective surface bright- 
ness for exponentiaUy corrected magnitudes as a function of 
bulge-t o-to tal ratio. The galaxies are the same type as described 
in Fig. The triangles represent the offset when seeing is not 
taken into account, and the squares show the effect of taking see- 
ing into account. There is a weak dependence on redshift. 



(A7) 



The probability of a galaxy of inclination lying be- 
tween i and i + di is: 



P{i < < i + di) = - sin(i)di 



(A8) 



Therefore the mean difference in measured and face-on 
effective surface brightness is: 



A /ie = 



Tr/2 



A^^o sm{i)di = -0.477 



(A9) 



Thus the calculated effective surface brightnesses of 
late-type galaxies may be 0.48 mag too bright. However 
more complicated effects such as seeing, the thickness of 
the disk and internal extinction will all tend to reduce the 
measured surface brightness of edge on disks to a greater 
degree than face-on disks, reducing the mean offset. 

The overall effects of bulge-disk decomposition and in- 
clination appear to be to make M* brighter by 0.1 mag and 
fil fainter by 0.55 mag arcsec"^, 0.5 mag arcsec"^ for Sas 
due to the bulge and 0.1 mag arcsec"^ for Sds due to the 
bulge and 0.5 mag arcsec"'^ due to inclination. It is difficult 
to be precise about this as the morphological mix of galaxies 
in the 2dFGRS is not known. 



where 

b — a cos(j) (A3) 

The surface brightness has increased at each point by 
a factor -^^^^^jj, increasing the semi-major axis until fj,{a) = 



a = Q0.41n(10)[At;im - A*o - 2.51ogjQ cos(j)] 

(A4) 

where a is the disk scale length. This increases the 
isophotal flux of the galaxy, as well as the isophotal radius. 



'ab 



(A5) 



mtot - 2.51ogig / 



(A6) 



An exponential profile is fitted to these parameters as 
in Cross et al. (2001). The central surface brightness and to- 
tal magnitude are calculated for this galaxy assuming that 
the galaxy is face on. The error in the central surface sur- 
face brightness is the difference between the true central 
surface brightness and the calculated central surface bright- 
ness. The error in the effective surface brightness is exactly 
the same, as the difference between central and effective sur- 
face brightness is a constant for an exponential profile. The 
total magnitude calculated above is the same as the true 
total magnitude. 
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